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Abstract 



A variant of a gauge theory is formulated to describe disclinations on Riemannian 
surfaces that may change both the Gaussian (intrinsic) and mean (extrinsic) curvatures, 
Q^ ' which implies that both internal strains and a location of the surface in R^ may vary. 

Besides, originally distributed disclinations are taken into account. For the flat surface, an 
^ ] extended variant of the Edelen-Kadic gauge theory is obtained. Within the linear scheme 

OO ] our model recovers the von Karman equations for membranes, with a disclination-induced 

^— ^ source being generated by gauge fields. For a single disclination on an arbitrary elastic 

^ . surface a covariant generalization of the von Karman equations is derived. 

o 

o\ 

1 Introduction 

d . Elastic two-dimensional structures which are free to change their geometry (membranes, thin 
^ i films, etc.) as well as deformable materials with spherical or tubular shapes (fullerenes, nan- 
otubes) are of considerable current interest (see, e.g., [0, ||, |^ and the references therein). The 
properties of these crystalline structures are found to be essentially affected by their topology. 
It has been found that topological defects, first of all disclinations, play an important role in 
these objects. In particular, the Kosterlitz-Thouless disclination unbinding transition in hex- 
atic membranes was shown to depend on shape fluctuations ^ since a membrane with a 
single disclination can lower its energy by buckling. Fullerenes and nanotubes always contain 
at least 12 disclinations on their surfaces (i.e. 5-fold coordinated sites) as a consequence of 
the Euler theorem. The effect of shape fluctuations on the interaction of the disclinations on a 
spherical surface with genus zero was studied in 0. 

Elastic models for membranes and shells are well known (see, e.g., [|^, H and the references 
therein). The main problem is how to incorporate defects into the elastic theory of two- 
dimensional fluctuating surfaces. A possible way has been considered in |l|] where defects were 
introduced ad hoc as source terms in the right-hand side of the von Karman equations. This 
approach is efficient in description of monolayers as well as membranes under condition that 
the bending rigidity which controls out-of-plane fluctuations is small. Similar approaches were 
developed in where the Coulomb-like model was formulated. Disclinations were introduced 
there as the point-like charges on the curved surface. 



1 



The modem trends in theoretical description of topological defects in condensed matter 



include geometrical and gauge-theory methods (see, e.g., 0, |T^, |T3l). In the present paper we 
put forward a gauge theory that enables one to describe disclinations on 2D elastic surfaces 
that may change both their intrinsic and extrinsic geometry. Namely, both the internal strain 
and the location of the surface in may vary. By analogy with the Edelen-Kadic (EK) gauge 
model of dislocation and dischnations 0, |15|, defects are incorporated via dynamical gauge 
fields. However, we formulate some principal statements which provide quite a new geometrical 
setting to handle a problem of describing defect-developing deformations of an elastic body. 
As is shown below, this allows to describe defects on arbitrarily fluctuating surfaces as well as 
dynamics of originally distributed defects. Notice that even in the 2D planar case our theory 
does not identically coincides with that obtained within the EK approach [1^. As a matter of 
fact, it includes the latter as a particular case and allows for a possibility to include objects 
with originally distributed disclinations. 

A basic motivation for our investigation has been that of extending the EK gauge theory 
to include nontrivial geometry in the 2D case, e.g., that of membranes, spheres, etc. It should 
be emphasized that within the standard EK approach there is no room for such objects at 
all. Indeed, within that approach deformation is mathematically described as a diffeomorphic 
mapping of a defect-free domain Dq G which is usually called a reference conflguration to 
another one D G -R" (current conflguration) with noneucledian metric tensor g. By the very 
deflnition Dq and D necessarily have the same dimension n, thereby ruling out the possibility 
to consider dynamics of, say, defects on 2D-surfaces for they are not domains in R^. 

Any attempt at generalizing the EK theory to study disclinations on general manifolds 
should imply two steps. First, one has to appropriately reformulate a classical theory of elas- 
ticity and, second, to introduce in thus obtained new geometrical setting dynamical gauge 
flelds. A straightforward generalization might be that of considering a diffeomorphic map 
X '■ Dq E —>■ D E Ti, where S stands for a Riemannian surface Q (i.e., a 2D real Riemannian 
manifold) located in R^. This map can (locally) be described by functions x°'{^'') {o-i^ = 1)2), 
where x"' denotes a point in Dq and corresponds to coordinates m D. A resulting Lagrangian 
to describe elastic properties of the media would follow as a function of the state vector x'^{.^)i 
which would result in a kind of a fleld theory in the nontrivial geometrical background. To 
include defects an appropriate gauge fleld on a curved space S would be further required, an 
obvious complication in contrast with the EK approach. 

The present approach has been proven to adequately incorporate dynamics of disclinations 
on varying elastic surfaces as well as to take into account originally distributed defects. The 
key observation that considerably simplifles and at the same time generalize the matter is that 
one should consider a surfaces S being embedded into a three-dimensional flat space R^ instead 
of proceeding entirely in terms of its intrinsic geometry. As we shall see shortly, following this 
idea will lead us to a fairly plain and surprisingly complete theory of defect dynamics on curved 
surfaces. 

The paper is structured as follows. In section 2 we formulate the gauge model of discli- 
nations on 2D elastic surfaces. The action which includes elastic deformations, self-energy of 
disclinations, and the curvature energy is constructed in a self-consistent way. A complete set 
of equations of motion is presented in section 3. To illustrate the model, we consider three ex- 
amples in section 4. First, we derive equations of motion for the planar case and show that they 
are distinct from those obtained within the EK approach. Second, we study a problem of the 



""^In particular, E may denote a Riemann surface, that is a ID compact orientable complex manifold. 
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fluctuating surface by employing the linear approximation. In this case, the known equations 
for fluctuating membranes are recovered in a self-consistent way with a source formed by the 
gauge disclination fields. Finally, a concrete realization of the model for a single disclination 
on arbitrary elastic surface is presented. Section 5 is devoted to concluding comments. 



2 The model 

Before proceeding, a few comments on the limitations of the theory as well as on the conventions 
employed are to be made. First of all, we are solely concerned with the 2D case, although our 
formalism can easily be extended to any space dimensions. Our motivation is that the 2D 
case is both the most important in applications and at the same time the simplest one in 
notation.0 Secondly, we take into consideration only the rotational symmetry of a system 
thereby making an attempt at describing disclinations and leaving aside dislocations and other 
defects. By turning to the full-fledged internal symmetry group instead of the orthogonal one 
the above restriction might be avoided, though it will evidently entail considerable conventional 
complications. Finally, static configurations are only considered, which does not seem, however, 
to be a deficiency in the following exposition, since including time evolution brings in no novel 
features compared to the EK theory. 

An action that is assumed to properly describe djTiamics of disclinations on a deformable 
elastic surface is taken in the form 



S — Sel + S gauge ~\~ Sfl, (1) 

where Sei describes the elastic properties of the media, Sgauge stands for the action of a gauge 
field that incorporates self-action for disclinations, and 5*/^ is the Helfrich-Canham action ||^, ^ 
to describe the energy of a free fluctuating surface. 

Let us start by discussing the first piece of the action. Let x°(a = 1,2) be a set of local 
coordinates on a certain Riemannian surface Eq. (Indices a, 6, c, ... = 1,2 are tangent to Sq, 
whereas i,j,k,... = 1,2,3 run over the basis of R^). Under a deformation. Eg is assumed to 
evolve into some other surface E. To describe this we find it convenient to introduce embeddings 
Eg, E — s> i?^ that can be realized in terms of two i?^-valued functions RIq^x^jx"^) and R^{x^ 



respectively. As the point is varied vectors i?(o) and R sweep surfaces Eq and E, 

respectively. This is nothing but a familiar two-parametric representation of surfaces in i?^, the 
point, however, being that 

%):=ri?(o) = i?(o)[0(a;)], (2) 

where 0* is a puUback of : Eq ^ E. In what follows functions i?(Q-)(x^, x^) are chosen to 
specify an initial configuration Eq, whereas dynamical variables R^{x^,x'^) are taken to describe 
the deformation Eq — E. With these conventions at hand, a proper generalization of the 
elasticity theory turns out to be a straightforward matter. 
Representations for the induced metrics follow immediately 

Qab = {g-Eo)ab = daR(0) ■ dbR(0), 



^11? case is irrelevant for us here since there are no II? objects with an intrinsic curvature 
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where the set {(p"'} stands for local coordinates on S. The strain tensor is determined to 



be \16, 17 



Eab — gab — gab- 

The elastic properties of the deformed surface are described by an action 

Sei = -\ dxHx^^{\{trEf + 2^xtrE^] , (4) 

where trE = g'^'^Eab, g = det\\gab\ \ and summation over repeated indices is assumed. Let us 
mention that we omit in (^) the terms of order E^ and higher. In some special cases they can 
be considered as well (see, e.g., [p!5[] ). 

To proceed, a properly defined gauge field to describe disclinations is to be introduced. To 
make this point transparent, let us step a bit aside and discuss the principle of a local gauge 
invariance from the geometric viewpoint. In this regard, a simple example may be of some help. 
Consider a scalar complex field '■ x G C. Let a Lagrangian exhibit a global U{1) 

symmetry: L ^ L under a transformation if) — > e*"?/', if) e~^°'ip. In what way the local gauge 
fields can then be introduced? To this end, one considers ip and ip as a description not of a 
map -R" —>■ C but of sections of a line C-bundle over i?" with the structure group U{1). This is 
a trivial bundle which admits global sections. A connection on this bundle ^ is a familiar gauge 
field Afj, which can be regarded as a U{1) valued one- form. The resulting theory is nothing but 
scalar electrodynamics. 

In order to incorporate disclinations in the elasticity theory (^), one should as the above 
example suggests consider the _R^-vector bundle over S and the i?^Q-)-bundle over Eq with the 

same structure groups S0{3). The so(3) valued one form A^^\x)dx°' {A'^^ = Wj^^ ■ L, L* G 
so(3)) serves as a connection one-form in the i?^Q-)-bundle space over Sq, with Wj^^^ being the 
gauge potentials. A connection on the i?^-bundle over Eq is obtained by pulling back the 
connection of the i?^-bundle over S: 

Wa:=<l>*Wak=da<l)\Wbk). 

By replacing in (||) ordinary derivatives daR and daR(o) by the covariant ones Va-R = daR + 

[Wa, R] and Va-R(o) = daR(o) + [W^i^'') -^{o)], respectively, one arrives at the desired locally 50(3) 
invariant representation for the elasticity Lagrangian. 

A few remarks are in order at this stage. First, we consider potentials Wj^^ as given fixed 
functions to describe disclinations originally distributed on Sq. These potentials being involved 
in a full set of equations of motion, provide a possibility to keep track of the dynamics of these 
disclinations under deformation. Such a possibility is missing in the standard EK theory where 
from the very beginning only defect-free initial configurations are allowed. 

Second, although essentially two-dimensional manifolds are considered, t/iree-dimensional 
rotational group S0{3) is involved, which seems to be quite natural in the framework of our 
approach. Enlarging the structure group to the semidirect product S0{3) > T(3), where T(3) 
stands for the group of translations in R^, enables one to include into consideration both 
disclinations and dislocations. 

Third, if we made an attempt to formulate an appropriate gauge theory in the scope of the 
direct generalization of the EK approach, we would encounter the following problem. With 

^To be more precise, a connection one-form is defined on the associated principal bundle P(i?"; U{1)) and 
completely specifies the covariant derivative on a C-bundle over i?". 
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the being not the S- valued functions defined on E but local sections of a S-bundle 

over S, the derivatives daXb must be replaced by suitable covariant derivatives. What would 
be the "gauge group" in the definition of these covariant derivatives? In general this would 
be the infinite dimensional group Diff E of all diffeomorphisms of S which is a fibre of the 
bundle. For an arbitrary S it almost surely does not include rotations in contrast to the group 
Diff = GL{3,R) 3 S0{3). To find a way out, one might be tempted to relate defects to 
a group of all isometrics of (S, (7) rather than to 5*0(3). This group belongs to Diff E and is 
generated by operators Va obeying 

^Va ^ = 0, 

where Cy^ denote Lie derivatives along corresponding vector fields Va- Needless to say that 
these equations cannot be in general solved explicitly, except in a few simple instances, e.g., 
planes, spheres and other highly symmetrical objects. In a general case, however, and as long 
as the S0{3) group is taken to be relevant for describing disclinations, this approach seems to 
pose a problem. 

Let us now turn to the second piece of the whole action, S gauge, which describes a self-energy 
of disclinations. It acquires a standard form of the 5*0(3) Yang-Mills action: 

Sgauge — T / ^s/ddx dx {J' , J' ab) ■, (5) 
4 J So 

where s is a coupling strength, the form (,) stands for the so(3) Killing trace and the so{3)- 
valued curvature tensor J^ab = Fab ■ L, Fab = daWb — dbWa + [Wa, Wb]- The so(3) generators 
Li obey the following commutation rules 

[Lj, Lj] = Lfc, j, = 1, 2, 3, 

where ef^ stands for the fully antisymmetric tensor in R^. 

Finally, the Helfrich-Canham action that describes a self-energy of a fluctuating two-surface 
in R^ looks like [g g 

Sfi = ^ f ^dx^dx^ {trKf + ^ / ^dxUx^ det g^^Kbe (6) 

where k is a bare bending rigidity and is a Gaussian rigidity, 

Kab = N- DaDbR (7) 
is the curvature tensor, and N is the unit normal to the surface 

fT_ [diR,d2R] 
\[diR,d2R]\' 

The covariant derivative 

Da:=da + Ta 

includes the Levi-Civita connection to be written down explicitly in the next section. Two 
scalar functions enter (^: trK = g°-^Kab called the mean (extrinsic) curvature, and 5 = 
detg"-^Kbe referred to as the Gaussian (intrinsic) curvature. In view of the fact that the second 
piece on the r.h.s. of (^ is a topological invariant and depends only on the genus of the surface, 
it does not affect classical equations of motion. Incorporating disclinations amounts then to 
performing in the above formulas a substitution 

da Va, 

which results in the observation that the both pieces of action (^) will now contribute to the 
equations of motion. 
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3 Equations of motion 

Equations of motion follow from the Hamilton principle of stationary action 6S = and read 

Vba'' + J = 0, (8) 



where a set of the stress vectors 
with 



sVaF'''-^[R,a'] + ^P' = (9) 

(T^ = i(V,i?)p»^ (10) 
p^^ := Xg^HrE + 2fiE''^ (11) 



has been introduced. The total covariant derivative 

Va := Va + 

includes the Levi-Civita (torsion-free, metric compatible) connection 

to take care of a metric factor ^ in Eqs. (^-|^) when a variation of the total action is calculated. 
It is noteworthy that depends on Wj^^ that enters in a definition of gab- 

+ (i^i°)i^r)4)-(w^i°^^(o))(H>ri?(o)). (13) 

We have also abbreviated 

f^J_^Sfi_ ^ ^ 1 SSfi 
SR ' VaSWb' 
bearing in mind that explicit formulae are available in particular cases. 

Obviously, for the elastic plane without defects (^ reduces to the well-known equilibrium 
equation dta^ = while (^ is absent. In the general case, both the gauge fields and the affine 
connection enter (||) and (H) thus affecting stress fields. In view of ([l2|) 

1 d 

and one may consequently rewrite Eqs. (§^ in the form 

1 

and 



^ a(v^?') + l»'i.,Sl + ./ = (14) 



^aJv^F-V^P'.fllH-^/'^O. (15) 

respectively. As is seen, the basic self-consistent equations of the model are strongly nonlinear 
and it is difficult to get the general solution of the problem. It will be shown below, however, 
that for some physically interesting problems these equations become essentially simpler and 
can be solved explicitly. 

^As was already mentioned, index a corresponds to the tangent space TEq, whereas i is referred to the basis 
of the underlying space R^. These are different manifolds and consequently a" = {af} are viewed as a set of 
vectors in R^. On the other hand, these spaces coincide in the EK approach and af reduces to a stress tensor. 
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4 Applications 



In this section we intend to consider three exphcit reahzations of the above theory. First, we 
show how the standard 2D EK theory emerges in the framework of our approach. Second, 
we consider an important in apphcations case of fluctuating elastic membranes, and, finally, a 
single disclination on arbitrary elastic surface is examined. 



4.1 Disclinations on elastic plane 

Let us start by recalling an explicit formulation of the 2D EK theory |T9|. Consider a diffeo- 
morphic map x '■ Dq ^ D, where Dq G appears as a strain free (which means gab = Sab) 
domain and D & stands for its image under deformation. Evidently, x° can be viewed upon 
as i?^-vector valued functions. Resulting action that describes dynamics of disclinations on a 
planar elastic body consists of two pieces Sei and S gauge given by Eqs. (|) and (j^), respectively, 
provided one puts 

Eab = Ba- Bb- Sab, Ba = daX + [Wa, x], 

where Wa = {0,0,Wa)- A complete set of the Euler-Lagrange equations can be easily derived 
and shown to possess exact vortex-like solutions |^ . 

On the other hand, basic notation of our theory in this case looks as follows. First we 

— * — * 

have the embeddings: i?(o) = {x,y,0) and R = {Ri{x,y), R2{x,y),Q), where x,y are Cortesian 
coordinates. Evidently, one has Wa = (0, 0, Wa), F"^ = (0, 0, F"''), where Fab = daWb - dbWa- 
The same representation holds for Wj^^^ which appears as a fixed gauge field associated with 
originally distributed disclinations. 

Equations of motion (H) and (|) take the form (J = = 0) 

-^db{^a') + [Wb,a'']=0 

l^da{y^F^') + ^Ja\R] = 0, 
where g = det\\gab\ \ and in view of Eq. ( |T3|) 

gab = Sab + eaaW^fi^ro) + ^^M'^ ^(0) + ( ) 4) , «,/? = !, 2. 

Metric tensor gab is seen to deviate from its fiat counterpart Sab at W^^^ ^ 0, whereby a 
nontrivial geometry is dynamically generated. It is at this point that our theory deviates from 
the standard EK one even in the trivial planar case. 

At W^^^ = the above equations reduce to those of the 2D EK theory In the linear 
approximation |2T| no interaction between Wa and Wj^'^^ fields occurs, so that a total disclination 



flow f yVdl that determines a source which gives rise to disclination-induced displacements, is 
found to heyV = W-W^^\ 



4.2 Disclinations on membranes 

Let us now turn to the gauge theory of defects for fluctuating membranes starting from (||,|D, our 
aim being flrst that to recover a conventional Landau theory without defects. No disclinations 
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are assumed to be originally distributed as well, so that we put WjP^ =0. A flat membrane 
that fluctuates in the 2;-direction can be described by the embeddings: 

^(0) = (a;,?/,0), R= {x + u^,,y + UyJ), 

where R = R^^^ + U and U = (ux, Uy, f) is a displacement of the (x, y, 0) point under deforma- 
tion. (In writing out the U components we stick to the Landau notation [^.) In accordance 
with the Landau theory U is assumed to be small in the transverse directions, along with the 
requirement that its z-component, f{x, y), is to be viewed as a slowly varying function. Precise 
restrictions on components of U directly follow those in It is clear by obvious reasoning 
that within the adopted approximations only the 2;-component of Wa matters, so that we may 
put Wa = (0, 0, Wa). The strain tensor is then determined to be 



Eab = Va-R ■ Vfe-R — 5ab 



with Greek indices a, /? = 1,2 being used to specify coordinates of the plane orthogonal to the 
2;-axes. 

As for the curvature tensor Kab, it is easily calculated to be 

Kab = dafdbf + 0{u\ udf, W^). (17) 
Within this accuracy one, consequently, obtains 

trK = A/. 

To proceed with the equations of motion, we write down Eq. (^ in components referred to 
the 2;-direction and transverse plane, respectively: 

VfeCTg + K J3 = 

and 

Vb(T^ + kJ^ = 0. 

Vector J is easily found to be J = (0, 0, A^/), which yields for the above equations 

^db{dafp''') + KA'f = 0, (18) 

and 

dbp'" = 0, (19) 



respectively. Equation (|Tl]) now reads 

p''^ = XtrES''^ + 2pE"'\ 



with E°'^ being given by (|T6|)- 

Within the present approximation, Wa = (0, 0, Wa), and Eq. @) is found to go over into 



daF'^' = ^^e.^pP'Rly (20) 
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In deriving this, we have also used that Jg = 0. Equations ( [I8|JT9| ) coincide exactly at Wa = 
with those in [0. In our case, however, there appears an additional equation ( pO|) which 
ensures the self-consistency of the model. As will be shown below, in the linear approximation 
this equation includes only gauge fields. In this case, one can use its solutions to determine the 
disclination-induced sources for the remaining equations. 

Notice that the foregoing procedure which allowed to obtain 
with the linearization scheme proposed within the EK model |0 



homogeneous scaling of the gauge group generators (see also details in [|T5 

1,^ one can get 



IS m agreement 
This scheme is based on a 
It can be shown 



that by applying this procedure to one can get (|18| -^0|). It is important to note, however, 
that in accordance with this procedure one has to choose properly the relation between the 
model parameters (A/s and fi/s). Depending on this choice one can describe elastic media with 
different properties. The classical elasticity theory which is of our interest here is recovered in 
the limit A/s ~ e and /i/s ~ e where e is a scaling parameter. In this case, (|20|) reduces to 



daF" 



0. 



A singular vortex-like solution of (|2l|) reads 

Wb = -u ebcdc log r, 

where u is the Frank index. 

One can easily see that ([l8| ) is rewritten as 



(21) 



(22) 



{dbdj)a 



ab 



(23) 



where (|T9|) is taken into account. Let us consider (0). It resembles the usual equihbrium 
condition but the stress tensor includes the gauge fields. Within the linear approximation 
terms with Wb can be separated thus forming the source in the right-hand side of (|19]). In 
particular, for the planar case one can reproduce the known exact solution for a straight wedge 



disclination 21 



To compare our results with those in let us differentiate Eq. (|l^), which yields 

dpdbp'" = 0. 

After straightforward calculations one can rewrite this equation as 

(A/4^ + 1/2) A trE = {d^dyff - dlfd^f + e^bdM. (24) 
Notice that the last term in the right-hand side of (0) describes a source due to disclination 



fields. For solution (UaI it takes the form 

eabdaWb = z/Alogr = 27rz/5(r). 

Introducing the Airy stress function x(^) by aba ='■ ^bm^andrndnXif) i oiie can finally rewrite 
equations (Bl) and (El) as 



(25) 



«:AV = {dlx){dlf) + {dlx){dlf) - 2{d,dyxmdyf), 
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respectively. Here Kq = 4/i(A + /i)/(A + 2/i), and trE = (1/(A + yu))Ax(r). As is seen, ( p5| ) 
are exactly the von Karman equations given in ^ for defects in hexatic membranes. It should 
be mentioned once more that the source term in ( |25D is not appeared ad hoc but generated by 
the gauge fields due to a disclination. Its exact form follows from the self-consistent solution 
of the basic model equations. An analysis of (^) shows that isolated positive (five-fold) 
disclinations on free membranes buckle into a cone, while the negative (seven-fold) disclination 
leads to a saddle surface. The energy of a positive disclination was found to be less than that of a 
negative one. It is interesting that this asymmetry is absent in flat membranes and monolayers. 
Notice that the linear approximation used in this section allows to properly describe only the 
small out-of-plane fluctuations. Otherwise, the full set of equations (Rpl) should be examined. 



4.3 Single disclination on arbitrary elastic surface 

In this subsection we demonstrate a non-trivial realization of the proposed model: a case of 
a single disclination on elastic surface. Let us consider a surface Eg in a three-dimensional 
Euclidean space R^. For any point p G Sq choose the 2;- axis to be normal to a tangent plane 
at p. Having in that way fixed the coordinate system in R^, we consider in what follows Sq as 
an embedding 

{u,v) {R'^Q^{u,v),RfQ^{u,v),Rlo^{u,v)), 

with = Ujx"^ = V being the local coordinates on Sq. Under a local deformation concentrated 
at p, any nearby point undergoes a displacement 

R(o){u,v) ^ R{u,v) = {Rfo-^{u,v) + U^{u,v), R\u,v)). 

As is above, the transverse displacement [/-*" is assumed to be small compared to the vertical 
one := R^ — -R^q)' latter being considered as a slowly varying function on a local chart 

that contains p. One may also put Wa = (0, 0, Wq), which is in agreement with the assumptions 
made. The strain tensor then becomes 

Eab = daRfo^dtU^ + e^pdaR^^,^Rfo)W, + {a^b) + daR^'d.W - daRl,^d,Rl,^ + 0{Ul, W\ V,_W). 

(26) 

In order to make a more close connection with the gauge theory of defects for fluctuating 
membranes of the preceding section, we find it appropriate to consider here a coordinate R^ as 
a dynamical variable rather than its displacement, f/^. In particular, the curvature tensor Kah 
appears then as a direct covariant generalization of Eq. (0). Indeed, it is clear that 

Nf:= [^--^(")'^--^(°)] =(0,0,1). 
|[f^«-R(o), <9t,i?(o)]| 1^ 

On the other hand, we have 

= iVp, ■ B,D^\, = {Np + 6Np) ■ DaD,Ri^, 

= Np ■ DaD,R\^, + O {\\5Np\\) , 6Np := Njy - N,. 

As Np are sufficiently smooth and slowly varying functions of a reference point p, there exists 
a certain neighborhood of p, Vp, such that ||5iVp|| ^ 1 for any p' G Vp. 
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In the linear approximation Np can be replaced by its unperturbed value N^^\ so that for 
a local deformation of Sq in the vicinity of p one obtains 

and consequently 

trK = g^'DaDbR' = DaD'^R' =: Aeo^,i^^ 

where the covariant Laplacian operator has been introduced. 

Equation of motion (|^) in the z and transverse components now reads 

1 



2 



{daRnp^' +KALi?^ = 0, (27) 



and 



(28) 



(0). 

respectively^. In the Monge gauge, i?(o) = (-R^o) = ='■ ^yR^o) = v =: y, R^Q^{x,y)), the last 
equation is written as 

Dbp"'' = 0. (29) 

As is seen, (^7D and (^9]) are nothing but a formally covariant representation of the Landau 
equations (|T8|) and ( pJ]) with the nontrivial geometry of Eg taken into account. These equations 
should be accompanied by the field equation (|^) which in the linear approximation describes 
the SO (2) gauge field in a curved background: 

DaF"'' = 0, F"'' = d^W^ - d^W. (30) 

The following steps are just the same as in the previous subsection. A singular solution of (|30| ) 
takes the form 

= -ve''^D,G{x,y), (31) 

where G{x,y) satisfies the equation 

DaD''G{x,y) = 27rS\x,y)/^, 

with Eab = \f9^ah being the fully antisymmetric tensor on Eq. Note also that a conventional 
flat" (5-function 5^(x, y) := 6[x)b{%)) is to be multiplied by a metric factor 1/y^ to make a 
scalar under a coordinate reparametrization. 
Equation (|27| ) now reads 

\{DbD^R^)p^' = -/^^ALi^^ (32) 
while (E9]) can be differentiated DaD^p"'^ = to yield 



)7 



(A/2/i + 1) Aeo. trE = {DaDbR'f - {/^co.R'f - {R' ^ i?(o)) + ^my6\x, y)/^, (33) 

where the condition e°'''DaWb = u/S.covG{x,y) has been taken into account. Obviously, (|33|) 
is the covariant analog of (^). It should be mentioned that due to the specific choice of the 
coordinate system related to a certain reference point, Eqs. (0) and ( p3D describe a single 
defect located at this point. In this regard, adding at least one more disclination makes the 
situation more difficult, so that turning to the basic equations of motion (^ and (^ seems 
necessary. 

^Components of the three- vectors i?(u, v) and R(q)(u, v) appear as scalar functions with respect to a change 
of coordinates on Sq and hence there is no need to replace the ordinary derivatives by the covariant ones in 
Eqs. (|7l H). 



11 



5 Conclusion 



The model developed in this paper allows to describe disclinations on arbitrary elastic sur- 
faces. It includes Riemannian surfaces that may change their geometry under deformations. 
In particular, within the proposed model one can study elastic properties of various materi- 
als containing disclinations: monolayers (flat surface), membranes (curved surface), fuUerenes 
(spherical surface) as well as nanotubes (which can be considered as deformed spheres). For the 
flat surface, we have obtained the extended variant of the EK gauge theory that includes orig- 
inally distributed disclinations. Within the linear scheme our model recovers the von Karman 
equations for membranes with a disclination-induced source being generated by gauge fields. 
For a single disclination on an arbitrary surface a covariant generalization of these equations is 
obtained. 

In our opinion, the most intriguing application of this theory might be that for a disclination 
on a sphere. The obtained equations (pTj) and (|29|) are the most general ones which allow to 
study this problem properly. Notice, however, that this is a challenging task. Indeed, an 



analytical solution of the simpler system (0) and (19) has been obtained only in the limiting 
case Kq —>■ oo. Nevertheless, there are precise numerical solutions of (|T^) and (|TU]) for arbitrary 
Ko (see, e.g., ||l|] and the references therein). Obviously, any attempts to solve ( pTj) and (|29|), 
either analytically or numerically would be of great interest. 

Another important problem relevant for the physics of fuUerenes and carbon nanotubes 
concerns the electronic properties of these materials. An attempt at extending the EK gauge 
theory of defects to include fermionic fields has been made in [^]. As has been shown, the 
self-consistent gauge model allows to describe physically interesting effects: the Aharonov- 
Bohm-like electron scattering due to disclinations [^, an electron localization near topological 
defects 0, as well as a formation of the polaron-type states near dislocations [0]. We 
expect that incorporating fermions in the above-formulated theory may provide a new insight 
into disclination theory in a curved background in the presence of electrons and reveal some 
novel physical phenomena. Fortunately, our approach allows for a natural extension of the 
model to include fermions. Indeed, the model action (|I|) is constructed out of sections of the 
vector bundles over Sq. Considering on the other hand a spin bundle over Sq, that is a tangent 
bundle TSq So with a structure group Spin{2) amounts to incorporating fermions into the 
theory, with fermionic fields being local sections of the spin bundle. (We assume that Sq admits 
a spin structure, which is the case, for instance, if Eq appears as a Riemann surface with genus 
g.) Since dimTjo = 2, the spin connection term drops out from an action, so that the Dirac 
operator on Eq takes the form 



where 7*^ = 6^7", gat = e"ef 6a/3 and the spin group Spin{2) is generated by two Dirac matrices 
7q, a = 1, 2 which can be taken as the Pauli matrices cti and a2- An explicit form for the fermion 
action as well as a full set of ensuing equations of motion will be given elsewhere. 

This work has been supported by the Russian Foundation for Basic Research under grant 
No. 97-02-16623. 
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